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SLIDING FREQUENCY STEADY- STATE PRECESSION MAGNETIC RESONANCE IMAGING 

Background of the Invention 

The present invention relates to the magnetic 
resonance imaging arts and particularly to steady- state free 
precession (SSFP) imaging. 

In steady- state free precession imaging, a series of 
magnetic resonance excitation pulses are applied at relatively 
short intervals, typically in the 2-15 millisecond range. 
After this series of pulses is applied for about three times 
the T x relaxation time (between one and six seconds for many 
common tissue types) , the magnetic resonance is stabilized in a 
steady-state condition. Thereafter, following each magnetic 
resonance excitation pulse, a phase -encoding gradient and a 
read gradient are applied to read out one phase -encoded data 
line until a full set of data lines are collected. 
Conventional steady- state free precession imaging sequences 
take advantage of signal generation along the longitudinal axis 
or both the longitudinal axis and the transverse plane. 

Area and moment balancing of gradient pulses which 
are applied in addition to frequency unwinding pulses which 
maintain coherence of the spin system rotating frame of 
reference have been used to increase image signal-to-noise 
ratios by making efficient use of the superposition of the 
rising and falling components of the resonance signal . 
Balanced gradient moments and short repetition times make 
steady- state free precession sequences insensitive to many 
types of patent motion and allow spins of flowing blood or 
cerebral fluid to come to a steady- state . 

Strong tissue contrast is obtained on the basis of 
the T x and T 2 relaxation times and the tip angle. Low tip angle 
images have a T x or proton density weighted contrast. High 
tip angle images have a contrast that is more T 2 dependent. At 
the Ernst angle, the image contrast is relatively independent 
of the repetition time and becomes proportional to the ratio of 
T a and T 2 of the imaged tissue. 
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Completely balanced sequences possess significant 
disadvantages. The pulse train produces a spectrally varying 
periodic excitation pattern with a frequency period 
proportional to the reciprocal of the sequence repetition time. 

This excitation pattern has two distinct null points per 
period that result in image banding artifacts. In order to 
achieve a stable steady- state in both the longitudinal and 
transverse magnetization, the excitation pulses are applied for 
several T x before the signal is sampled. Often, the 
stabilization time is a significant fraction of the total 
imaging time, often more than half. 

The periodic nature of the excitation pattern leads 
to banding artifacts that tend to follow the magnetic field 
contours in steady-state free precession images, because 
magnetic field inhomogeneity contributes to the local spectral 
offset in each imaging pixel. In general, the appearance of 
banding artifacts is governed by the total phase precession 
gained by spins during the interpulse (TR) period. Common 
sources of phase precession in liquid state imaging include: 
resonance offset, chemical shift, magnetic field inhomogeneity, 
spin coupling, unbalanced gradient pulses, and unbalanced 
changes in the transmitter carrier frequency and receiver 
demodulation frequency or their phases. An effective (average) 
spectral offset (radians/sec) can be defined for each spin 
species in any pixel by dividing its total phase gain between 
pulses (radians) , by the interpulse spacing (seconds) . 

One proposed solution for banding artifacts is the 
use of very short repetition times to reduce the total phase 
gain from local changes in spectral offset. Another solution 
is to average a plurality of excitations or images with a range 
of spectral offsets such that the banding artifacts average out 
or cancel, because of the periodic nature of the steady-state 
spectral response. Both solutions, however, have drawbacks. 

Ultra-short repetition times compact the time 
available for phase-encode and read gradients. As the times 
become shorter, the gradients need to turn on and off more 
quickly, i.e., a higher slew rate. Typically, higher strength 
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gradients are used to offset the shorter read times. Such 
systems are expensive and may be unfavorable to patient safety 
and comfort. Rapidly changing magnetic fields can induce both 
muscular and neural stimulation. The use of strong gradients 
5 also requires the use of wide bandwidth receivers to digitize 
the read out resonance signals. Such wide bandwidth receivers 
decrease the signal -to-noise ratio counteracting a primary 
advantage of the steady- state free precession method. Further, 
short repetition times severely restrict image resolution. 

10 Averaging spectral offset magnetic resonance signals 

or complex valued images created from such signals successfully 
cancels banding artifacts. The RF pulses in the steady-state 
free precession excitation train are shifted in constant phase 
increments, e.g., 90°, for each successive RF pulse. A similar 

is phase shift is also applied to the signal detected in the 
imagerOs receiver. These combined phase shifts change the 
frequency of the rotating reference frame. Hence, a 90° phase 
increment moves the spectral excitation pattern by a frequency 
offset equal to 1/4 of a period. A 180° phase increment shifts 

20 the pattern by 1/2 of a period, and so on. Collecting a set of 
data lines for each of a plurality of different phase cycles 
with different phase increments spaced equally over any 3 60° 
period allows averaging to cancel the banding artifacts. 
Alternately, a Fourier transform applied across the spatially 

25 encoded data dimension can be used to average or combine the 
phase shifts. The collection of the steady-state free 
precession signal with a series of increment phase cycles is 
the equivalent to encoding the data with spectral offset in 
addition to the usual spatial encodings that are accomplished 

30 by the gradient pulses. Calculation of a mean or averaged 
signal is equivalent to the calculation of the central (k=0) 
Fourier coefficient for the spectral domain. 

A primary drawback of the signal averaging or 
transform approach is that a large number of averages are 

35 collected to cancel the banding artifacts, typically at least 
four. Time is allocated for the system to regain its steady - 
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state after each change of the phase cycle. This causes image 
data to be averaged in a serial fashion rather than in a 
sequential fashion. That is the four or other plural number of 
like data line sets with each of the averaged spectral offsets 
are collected at temporally different times rather than 
sequentially. 

The present invention provides a new and improved 
steady- state free precession imaging technique that overcomes 
the above -referenced problems and others. 

Summary of the Invention 

In accordance with the present invention, a method of 
magnetic resonance imaging includes placing a region of 
interest of a subject in an examination region in a main 
magnetic field. Resonance is repeatedly excited with an RF 
pulse at regular repeat intervals until a steady-state 
condition is established. After the steady- state condition has 
been established, RF pulses continue to be applied and at least 
phase-encode and read gradient pulses are applied following 
each RF pulse to acquire a set of data lines in k-space. With 
the data acquisition of each data line, a spectral offset is 
changed to modulate the steady-state condition sufficiently 
slowly that the steady- state condition does not need to be 
reestablished between the acquisition of consecutive data 
lines. Steady-state signals derived from a plurality of the 
data lines are combined. 

In accordance with another aspect of the present 
invention, a steady-state free-precessing magnetic resonance 
imaging sequence is provided. Magnetic resonance is excited to 
establish a steady-state condition and a plurality of sets of 
data lines are collected with different spectral offsets. The 
spectral offsets are changed sufficiently slowly that the 
steady- state condition does not need to be reestablished 
between the acquisition of consecutive data lines. 

In accordance with another aspect of the present 
invention, a magnetic resonance image generating apparatus is 
provided. A means repeatedly excites resonance in a region of 
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interest of a subject in an examination region in a main 
magnetic field with RF pulses at regular intervals until a 
steady-state condition is established. The means also 
continues to apply the RF pulses after the steady-state 
condition has been established. A means applies at least 
phase -encode and read gradient pulses following each RF pulse 
to acquire a set of data lines in k- space. A means changes a 
spectral offset with the generation of each data line 
sufficiently slowly that the steady-state condition does not 
need to be reestablished between the acquisition of consecutive 
data lines. A means combines the steady-state signals derived, 
from the plurality of data lines. 

One advantage of the present invention is that it 
reduces data acquisition time. 

Another advantage of the present invention is that it 
reduces or eliminates banding artifacts. 

Another advantage of the present invention is that it 
is amenable to higher resolution imaging. 

Still further advantages of the present invention 
will become apparent to those of ordinary skill in the art upon 
reading and understanding the following detailed description of 
the preferred embodiments. 

Brief Description of the Drawings 

The invention may take form in various components and 
arrangements of components, and in various steps and 
arrangements of steps. The drawings are only for purposes of 
illustrating a preferred embodiment and are not to be construed 
as limiting the invention. 

FIGURE 1 is a diagrammatic illustration of a magnetic 
resonance imaging system in accordance with the present 
invention; 

FIGURE 2 is a diagrammatic illustration of a segment 
of a steady- state free precession imaging sequence; 

FIGURE 3 illustrates timing of changes in the RF 
phase increment and the amplitude of the phase -encode gradient 
pulses; 
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FIGURE 4A illustrates timing of changes in the RF 
phase increment and the amplitude of the phase -encode gradient 
pulses with more data lines at the center of k- space being 
collected and averaged than adjacent edges of k- space; 

FIGURE 4B illustrates the regions of k- space that are 
collected in each cycle portion; 

FIGURE 4C illustrates the number of redundant data 
lines collected for various regions of k- space; and, 

FIGURE 5 illustrates yet another embodiment in which 
each set of data lines spans about half of k-space. 

Detailed Description of the Preferred Embodiments 

With reference to FIGURE 1, the present invention is 
equally applicable to high field and low field systems. It is 
15 also applicable to both open or C-type magnets and cylindrical 
bore type magnets, although an open magnet is shown in FIGURE 1 
for purposes of illustration. An imaging region 10 is defined 
within a B 0 field generated by superconducting or resistive 
coils 12a, 12b of a main field magnet. The main field magnet 
20 generates a temporally constant, main magnetic field B 0 through 
the imaging region along a vertical axis, in the open magnet 
embodiment. A subject is supported in the imaging region by a 
patient support 14. 

Gradient magnetic field windings or coils 20a, 20b of 
25 a gradient field coil set are disposed adjacent the imaging 
region 10 to create magnetic field gradients across the main 
magnetic field, along three orthogonal axes. The gradient 
coils are connected to gradient amplifiers 22 and a gradient 
magnetic field controller 24. In the illustrated open magnetic 

30 system, the gradient coils are mounted above and below the 
patient. In bore type systems, the gradient coils are 
traditionally mounted on a tube that surrounds the bore. 

In order to excite magnetic resonance in selected 
dipoles of the subject disposed in the imaging region, windings 
35 or coils 30a and 30b of an RF coil are disposed between the 
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subject and the gradient coils. Although whole body RF coils 
are illustrated, it is to be appreciated that head and other 
localized or specialized coils may also be utilized. At least 
one radio frequency transmitter 32, preferably a digital 
transmitter, causes the radio frequency coil to transmit radio 
frequency pulses with the amplitude, timing, and phase set by a 
radio frequency pulse controller 34 into the imaging region. 

With continuing reference to FIGURE 1 and further 
reference to FIGURE 2, a sequence controller 36 causes the 
gradient controller 24 and the RF controller 34 to implement a 
steady-state free precession technique. More specifically, a 
series of RF or radio frequency excitation pulses 40 are 
applied displaced by an interpulse time interval or repeat time 
TR. In the preferred embodiment, the repeat time is about 2- 
15 milliseconds. A train of RF pulses are applied to establish 
the steady state. Once the steady state is established, a 
phase -encode gradient 42 is applied to select a data line in k- 
space between each pair of RF pulses. A phase-encode rewind 
gradient 42' removes the phase -encoding after each data line is 
read. Also between each pair of pulses, a read gradient 44 is 
applied for frequency- encoding a resultant magnetic resonant 
echo 46. In one embodiment, a slice select pulse 48 is applied 
in conjunction with each RF pulse to limit the excitation to a 
preselected slice. 

By varying the frequency of the slice select RF pulse 
by a means consistent with maintaining the coherence of the 
spin system rotating frame then the physical location of the 
excited slice may be stepped to generate multiple two- 
dimensional images covering any desired volume. The slice 
select frequency shift may be created for a constant 
transmitter carrier frequency by modulating the stored 
specification of the selective RF pulse in the transmitter 
waveform memory with an audio frequency equivalent to the 
frequency shift of the slice from the magnet center when the 
slice select gradient is applied. Alternatively, the slice 
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select frequency shift may be created by shifting the 
transmitter carrier frequency during transmission of the RF 
pulse in a phase continuous fashion. The phase difference 
between transmitter and receiver channels must then be rewound 
when the transmitter is returned to the central . magnetic 
resonance frequency after transmission of the RF pulse. Phase 
rewinding between the transmitter and receiver can be 
accomplished by setting the transmitter frequency in the 
opposite sense to that used for slice selection for an 
equivalent time interval to the duration of the RF pulse. 

Alternately, a second phase-encode gradient can be 
applied along the ""slice select direction for three- 
dimensional data acquisition with two dimensions phase-encoded 
and one dimension frequency encoded. 

The magnetic resonance signals 46 are picked up by 

the whole body RF coil 30a, 30b or a localized coil (not 
shown), amplified, and demodulated by a receiver 50, preferably 
a digital receiver. The transmit and receive coils can be the 
same coil or different coils. 

Initially, the RF pulses 40 can be applied, either 

with or without the imaging gradients to excite magnetic 
resonance repeatedly and bring the spin system of the subject 
into a steady-state. Typically, a time period of about three 
T x is sufficient to establish a steady-state. For body tissue 
with a T x of about 2, 000 milliseconds, the RF pulses are 
applied for about 6,000 milliseconds, i.e., 6 seconds. With a 
repeat time TR of 6 milliseconds, 1,000 RF pulses are applied 
to establish steady-state. 

Once steady- state has been achieved, the phase-encode 
gradient 42 and read gradient . 44 are applied between each 

adjacent pair of RF pulses. 

In the illustrated embodiment, the phase-encode 
gradient 46 is incremented in each repetition along a 
trajectory 52 to step through k- space. Although the trajectory 
52 is shown as a continuum in FIGURE 3, it is to be understood 
that it represents the amplitude change of the 128, 256, etc., 
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phase-encode pulses that span k-space and not a continuous 
phase-encode gradient. When collecting an image with 128 lines 
in k-space, there are 128 repetitions of the RF pulse 40 
between a beginning 54 and an end 56 of the phase -encode 
gradient trajectory 52. For a repeat time of 6 milliseconds, 
128 data lines are collected in just under 800 milliseconds. 

Once the steady-state has been established, an RF 
phase adjusting means or processor 60 adjusts the phase of the 
RF transmitter and receiver in coordination. The phase of the 
transmitter and the receiver are incremented just prior to 
transmission of each RF pulse and held constant through the 
rest of the interpulse interval. Although the trajectory 62 is 
shown as a continuum in FIGURE 3, it is to be understood that 
it represents a series of phase increments that span a period 
of 360°. It is desirable to keep changes between successive, 
phase increment steps less than about 4° so that the steady- 
state condition of the magnetization can reversibly anneal 
whilst the effective spectral offset is modulated. 

The trajectory of the phase-encode gradient and a 
trajectory 62 of the RF phase are coordinated such that over 
the course of collecting one of N sets of data lines between 
points 54 and 56, the RF phase increment has shifted 3 60°/N, 
where N=4 in the embodiment illustrated in FIGURE 3. When a 
full set of data of k-space is collected, the phase-encode 
gradient magnitude is returned to the initial value and the 
process is repeated to generate a second set of data lines 
spanning k-space. 

In an illustrated embodiment in which averaging takes 
place in k-space, a data line sorting means or processor 70 
sorts the data lines among a series of memories 72 such that 
the first data line of each of the N sets is stored in a first 
memory 72 x , the second line of each of the N sets is stored in 
a second line data memory 72 2 , and so on. Once N corresponding 
data lines have been collected, a summing or averaging 
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processor 74 averages the like data lines which combined data 

lines are then loaded into a k-space memory 76. An inverse 

Fourier transform processor 78 Fourier transforms the k-space 

data to create an image representation which is stored in an 
5 image memory 80. Preferably, the image representation is a 

three-dimensional representation and can be generated either as 
a series of contiguous slices or by a three-dimensional Fourier 
transform. A video processor 82 accesses the image memory to 

generate slice, surface rendering, angiographic, or other 
10 images for display on a monitor 84 . 

Although the phase-encode magnitude is illustrated in 
FIGURE 3 as being linear, i.e., that the lines of k-space are 
collected serially, it is to be appreciated that these data 
lines can be collected in other orders. Moreover, some of the 

is lines of k-space can be collected more than once in each of the 
N segments of the imaging cycle. Further, the rate at which 
the spectral offset (RF phase increment) changes need not be 
linear. The change can plateau, or change more or less rapidly 
at various points in the data collection cycle. For example, 

20 the RF phase increment might change more quickly near the edges 
of k-space and more slowly during the collection of data lines 
near the center of k-space. 

The embodiment shown in FIGURE 3 cancels banding 
artifacts by combining the theoretical minimum number of 

25 complete data sets needed for this purpose (N=4) . Such a set 
of data can be said to form a complete artifact cancellation 
set. Better cancellation of banding artifacts is obtained if 
the size of the cancellation set N is increased, because this 
provides better sampling over the spectral response of the 

30 steady-state signal. While this is a desirable improvement, it 
comes at the expense of scan time. 

The cost of artifact cancellation can be reduced 
because it is unnecessary to combine complete data sets. Only 
subsets of the complete data near the center of k-space need be 

35 acquired if at least four full k-space memories have been 
filled to provide enough data for basic artifact cancellation 
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in every line of k-space. Improved artifact cancellation can 
be obtained with fewer data lines, because band artifacts tend 
to be most strongly represented in the central data lines of k- 
space which contribute most to the overall image contrast and 
intensity, while the outer data lines of k-space contribute 
mainly edge definition and in large part the image noise. 

One illustrated embodiment of the invention that 
takes advantage of these principles to provide a scan time 
saving equal to a factor of two is shown in FIGURE 4A. Once 
steady- state has been established, the RF phase of the 
transmitter and receiver devices is incremented in steadily 
increasing steps from 0° to 360° and four full k-space data sets 
are acquired in the same fashion as is shown in FIGURE 3. Once 
the RF phase increment reaches a value of 360°, a second period 
of RF phase incrementation can begin. While it might appear 
that there is an extreme discontinuity in the RF phase 
trajectory, in reality there is none, because of the 
equivalence of 0° and 360° phase shifts. 

During the second period of RF phase incrementation 
shown in FIGURE 4A, twelve data sets covering the central 
region of k-space are acquired. The states of the data 
memories 72 are shown in FIGURE 4B and their resultant 
combination after the averaging processor 74 is shown in FIGURE 
4C. The twelve central regions of k-space shown in FIGURE 4B 
can be grouped into three subsets. The order in which the data 
lines in each subset are acquired relative to the RF phase 
trajectory (spectral offset) is chosen so that each subset 
forms a complete cancellation set for the region of k-space 
acquired . 

For the embodiment shown in FIGURE 4A, the first 
subset contains three eights of the data lines needed for a 
complete k-space data set, two thirds of these data lines are 
acquired before the k=0 point and one third of them after. The 
second subset also contains three eights of the data lines 
needed for a complete k-space data set, one third of these data 
lines are acquired before the k=0 point and two thirds of them 
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after. The third subset contains the data lines in the central 
quarter of k~ space, FIGURE 4C shows the total number of data 
lines that contribute to various regions of k-space when all 
the data for the scan shown in FIGURE 4A is combined. The 
5 central quarter of the combined k-space results from sixteen 
data lines, the two flanking regions each spanning one eighth 
of combined k-space results from the combination of eight data 
lines, and the two outer quarters of k-space result from the 
combination of four data lines. This unevenness in the 

10 combined data lines can be compensated for by applying a 
filter, the simplest of example which normalizes each data line 
in the combination by the number of individual data lines from 
which it is derived. The overall result of the acquisition is 
an image which is substantially equivalent in quality to that 

15 obtained using N=16 with the embodiment shown in FIGURE 3, but 
acquired in half the time. It is to be appreciated that the 
number of data sets used in the cancellation subsets N, the 
size and extent of the central regions of k-space, their 
ordering with respect to the RF phase increment scheme, and the 

2 0 RF phase increment scheme used in the embodiment illustrated in 
FIGURE 4A are merely selectable properties of a particular 
scan, or scanning scheme, and that many different variations of 
this particular embodiment with different tradeoffs between 
scan time and scan quality are possible. 

25 With reference to FIGURE 5, in each phase-encode 

subcycle, the data lines spanning just over half of k-space are 
sampled. That is, the data lines near the center of k-space 
are sampled each time; but in alternate samplings, data is 
sampled from the central lines only to one side of k-space. 

30 Preferably, the data on each side of k-space are sampled 
alternately. In a second 360° of RF phase rotation, data 
collection is preferably offset by 360°/2N. Each portion (just 
over half) of the data lines of k-space is combined or averaged 
with the set of data lines from the opposite half of k-space 

35 collected at the closest RF phase angular offset, i.e., closest 
to 180° apart or closest to 0°, 360°, etc., apart. Appropriate 
filtering is used to correct for the redundancy of data at the 
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center of k-space. If each data line is averaged with the 
closest other half leading and trailing data line, appropriate 
filtering is applied for proper weighting. 

The half k-space data line collection scheme, like 
5 other collection schemes, can be formed at 90° intervals, 45° 
intervals, 60° intervals, or the like. For purposes of 
illustration, the embodiment of FIGURE 5 is illustrated with 
the first 360° cycle of RF phase being collected at 60° 
intervals, i.e, 0°, 60°, 120°, 180°, 240°, and 300° and the 

10 second 360° cycle of the RF phase is collected at 30° offset 60° 
intervals, i.e., 30°, 90°, 150°, 210°, 270°, and 330°. In this 
manner, data collected at 12 angular increments is averaged for 
more complete artifact cancellation. 

Rather than averaging the half k-space data sets 7 the 

15 remaining half can be completed by conjugate symmetry. 

It is to be appreciated that the partial acquisition 
techniques described above can be used with the conventional, 
multiple spectral offset method as described in Zur, et al . , 
""Motion Insensitive, Steady-State Free Precession Imaging^ , 

20 Mag. Res. Med. 16:444-459 (1990) rather than the preferred 
continuous spectral offset shifting. Further, the RF phase 
increment can change in a piecewise linear manner or a 
continuously varying manner. For example, the RF phase 
increment can change along a continuum adjacent the edges of k- 

25 space and plateau and remain constant for data lines taken near 
the center of k-space for improved matching and band 
cancellation. This piecewise linear RF phase increment change 
still eliminates the restabilization periods that would occur 
between stepped RF phase increments needed to acquire data at 

30 multiple spectral offsets. 

Although this technique has been described with 
reference to two-dimensional imaging for simplicity of 
illustration, it is to be appreciated that it is equally 
applicable to three-dimensional imaging techniques. 

35 In addition to slow modulation of the spatial offset 

frequency, the RF pulsing may also have a modulated irradiation 
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with components at multiple spectral offsets, i.e., a 
"chirped 1 1 RF pulse train. In one embodiment, multiple trains 
of identical RF pulses are cycled each with its own phase shift 
incrementation scheme being applied to the subject in a 
synchronous fashion. As the RF pulses occur at exactly the 
same time, they can be combined as complex phasors to obtain a 
single RF pulse with a new amplitude and phase at each position 
in the steady- state free precession train. In this manner, 
more than one spectral frequency can be simultaneously 
irradiated. As the amplitude of the RF pulse train components 
are varied, there is effectively a superposition of complex 
excitation trains. A very large number of sequences can be 
generated from RF pulse trains formed by summing complex 
contributions of the form: 

Ae i(am+P) 

where A is the amplitude that controls the relative strength of 
irradiation at a particular spectral frequency, i is the square 
root of -1, co is the phase increment used to cycle the n th pulse 
measured in radians. The parameter co sets the spectral offset 
of the pulse train relative to the main carrier. The parameter 
n is the index of the n th RF pulse in excitation trains 

n=l,2,3, The parameter p is an overall phase shift applied 

to a component excitation train where 0^p<27c radians. Some 
sequences generated by this method can be reduced by the 
removal of unneeded slice selection gradients to being the 
equivalent of a basic steady-state free precession sequence 
with multiple sampling windows placed in the precession period. 

In some manifestations, some of the RF component amplitudes 
become zero . 

A useful class of sequences can be obtained by 
restricting the amplitude factors A of the input trains to be 
equal, choosing the phase increments co to generate equally 
spaced bands at higher multiples of the periodicity set by the 
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repeat time TR, in selecting the overall phase of each 
component of radiation so that the final combined amplitude 
numbers of each RF pulse in the train are equal . 

The spectral offset modulation scheme generated by- 
such a combination of conditions causes the instantaneous 
spectral offset with which each RF pulse is transmitted to vary 
cyclically through a finite set of values and creates a class 
of sequences that generates a periodically repeating cycle of 
steady-state signals, each of which has a different frequency- 
shifted spectral response pattern. During imaging, the signals 
generated for each cycle are each encoded with the same phase 
encode state and summed together. It will be appreciated that 
this converts the imaging operation from a serial averaging 
process into a sequential one. 

The spectral responses formed in the multiple 
irradiation case are periodic in offset frequency, having a 
banding pattern in longitudinal and transverse magnetization. 
Conventional steady- state free precession has one irradiation 
frequency and generates two intense bands per cycle of spectral 
offset in its periodic response. Irradiation with two offset 
frequencies generates summed patterns with four less intense 
bands per cycle. Irradiation with three, four, etc. 
frequencies generates six, eight, etc. bands per cycle, 
respectively. The banding artifacts possess the characteristic 
of having two spectral periods, a major one (low spectral 
frequency) set by the choice of interpulse spacing, and a minor 
one (high spectral frequency) set by the number of irradiation 
frequencies in the RF pulse train. 

With conventional balanced SSFP imaging, the 
interpulse spacing or repeat time TR is typically made shorter 
to minimize banding artifacts. However, when multiple spectral 
offsets are irradiated, increasing the repeat time TR 
facilitates the condition at which the bands due to the minor 
spectral response period are pushed to a size of less than one 
pixel and the band and non-band information within each pixel 
is averaged to provide a uniform image without banding. The 
effect of signal variation due to the major period of the spin 
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system spectral response ensures that the local average over 
any minor period is never zero. 

One feature of the technique with a large number of 
irradiation components is it exhibits an anomalous time-of- 
flight effect and that slow moving blood appears bright while 
fast moving blood appears black in imaging. Stated in more 
practical terms, the images have bright veins and dark arteries 
while the stationary tissue is an intermediary gray scale. 

Better cancellation of band artifacts may be achieved 
with these sequences using a relatively small number of 
averages. In part, band artifacts are not as strong when 
multiple spectral offset irradiation is used. Also in part, 
the frequency of banding increases with the number of 
independent irradiation sequences. If the bands are so closely 
spaced that they fall within the natural resonance offset 
distribution within each imaging pixel, then the steady- state 
free precession banding artifacts disappear automatically. 
This effect is controlled by the choice of the repeat time. 
Increasing the repeat time TR will lower the band frequency 
spacing and lead to internal averaging within pixels. For four 
irradiation frequencies and a 5 millisecond repeat time, the 
bands are spaced 50 hertz apart, which is close to the value 
for obtaining internal pixel averaging in humans. The total 
imaging time for this sequence would be identical to that to 
obtain four averages with the conventional sequence. 
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Having thus described the preferred embodiment, the 
invention is now claimed to be: 

1. A method of magnetic resonance imaging 
comprising: 

placing a region of interest of a subject to be 
imaged in an examination region in a main magnetic field; 

repeatedly exciting resonance with an RF pulse at 
regular repeat intervals until a steady-state condition is 
established; 

after the steady-state condition of a magnetic 
resonance signal has been established, continuing to apply the 
RF pulses and applying at least phase-encode and read gradient 
pulses following each RF pulse to acquire a set of data lines 
in k- space; 

with the acquisition of each data line, changing a 
spectral offset to modulate the steady-state condition 
sufficiently slowly that the steady- state condition does not 
need to be reestablished between the acquisition of consecutive 
data lines; 

combining steady-state signals derived from a 
plurality of the data lines. 

2. The method as set forth in claim 1, wherein the 
spectral offset is changed sufficiently slowly that the steady- 
state condition is reversibly modulated between the acquisition 
of each data line, 

3 . The method as set forth in either one of claims 
1 and 2, wherein: 

the spectral offset is changed sufficiently slowly 
that the steady- state condition is not destroyed. 

4. The method as set forth in any one of claims 1- 
3, wherein, in the repeatedly exciting step, the RF pulse train 
has at least two spectral offset components. 
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5. The method as set forth in any one of claims 1- 
4; wherein the time integral of the magnetic gradients is zero 
over the repeat interval . 

6. The method as set forth in any one of claims 1- 
5, further including: 

combining corresponding data lines of a plurality of 
sets of data lines in k-space collected with different spectral 
offsets; 

reconstructing the combined corresponding data lines 
into an image representation. 

7. The method as set forth in any one of claims 1- 
5, further including: 

reconstructing each set of data lines into an image 
representation ; 

combining image representations generated with data 
lines collected with different spectral offsets. 

8. The method as set forth in any one of claims 1- 
7, wherein data lines near a center of k-space are sampled more 
often than data lines near an edge of k-space, over a plurality 
of sets of data lines. 

9 . The method as set forth in claim 8 wherein each 
set of data lines covers a central region of k-space and half 
of a remainder of k-space. 

10. The method as set forth in claim 9 further 

includes : 

generating data lines for the remaining half of k- 
space from complex conjugates of the data from the collected 
half of k-space. 

11. The method as set forth in either one of 
claims 9 and 10 , further including: 
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combining sets of data lines which span one half of 
k-space with sets of data lines which span the other half of k- 
space but collected with different spectral offsets. 

12. The method as set forth in any one of claims I'- 
ll, wherein in the repeatedly exciting step an RF phase 
increment is applied to successive RF pulses and varied over a 
range of 360° while N sets of data lines are collected, each 
set of data lines being collected with RF phase increments 
offset by 3 60°/N relative to adjoining sets of data lines. 

13 . The method as set forth in claim 12 wherein the 
RF phase increment changes continuously and linearly from 
commencement of collecting data lines of a first set of data 
lines after establishing steady- state free precession through 
collection of a last data line of a last set of data lines 
collected. 

14. The method as set forth in claim 12 wherein RF 
pulses are applied to establish steady-state free precession 
over a duration at least the shortest nuclear relaxation time 
of material in the imaged region of interest of the subject. 

15. The method as set forth in any one of claims 12- 
14, wherein a full data line is collected between adjacent RF 
pulses and each data set includes at least four data lines. 

16. A steady-state free precession magnetic 
resonance imaging sequence in which magnetic resonance is 
excited to establish a steady- state condition and a plurality 
of sets of data lines are collected with different spectral 
offsets, characterized by: 

changing the spectral offset sufficient slowly that 
the steady-state condition does not need to be re-established 
between the acquisition of consecutive data liens. 
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17. The steady-state free-precession magnetic 
resonance imaging sequence as set forth in claim 15 wherein the 
spectral offset of the sequence sufficiently slowly that the 
steady-state condition is reversibly modulated between the 
acquisition of each data line. 

18. The steady-state free precession magnetic 
resonance imaging sequence as set forth in either one of claims 
16 and 17, further characterized by: 

changing the spectral offset continuously. 

19. The steady-state free precession magnetic 
resonance imaging sequence as set forth in either one of claims 
16 and 17, further characterized by: 

changing the spectral offset piecewise continuously. 

20. The steady-state free precession magnetic 
resonance imaging sequence as set forth in either one of claims 
16 and 17, further characterized by: 

changing the spectral offset piecewise linearly. 

21. The steady-state free precession magnetic 
resonance imaging sequence as set forth in any one of preceding 
claims 16-2 0, wherein changing the spectral offsets includes: 

shifting RF phase. 

22. The steady-state free precession magnetic 
resonance imaging sequence as set forth in any one of preceding 
claims 16-20, wherein changing the spectral offset includes: 

shifting the RF frequency. 

23. The steady-state free precession magnetic 
resonance imaging sequence as set forth in any one of preceding 
claims 16-20, wherein changing the spectral offsets includes: 

sweeping the magnetic field. 
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24. The steady-state free precession magnetic 
resonance imaging sequence as set forth in any one of preceding 
claims 16-23, further characterized by: 

exciting magnetic resonance with a radio frequency 
pulse train with at least three spectral offset components. 

25. The steady-state free precession magnetic 
resonance imaging sequence as set forth in claim 24, further 
characterized by: 

combining the periodically modulated steady- state 
signals generated by an integral number of periods of the RF 
excitation pulse train. 

26. A steady-state free precession magnetic 
resonance imaging method comprising: 

exciting magnetic resonance concurrently at a 
plurality of frequencies; 

establishing a periodically modulated steady-state 
condition of magnetic resonance signals; 

combining the periodically modulated steady- state 
signals generated by an integral number of periods of the RF 
excitation pulse train; 

repeatedly sampling the combined resonance signal to 
generate a plurality of data lines; 

reconstructing the data lines into an image 
r epr e s ent a t i on . 

27. An apparatus for generating magnetic resonance 
images comprising: 

a means (32 , 34) for repeatedly exciting resonance in 

a region of interest of a subject in an examination region in a 
main magnetic field with RF pulses at regular repeat intervals 
until a steady-state condition is established and for 
continuing to apply the RF pulses after the steady-state 
condition of a magnetic resonance signal has been established; 
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a means (22 , 24) for applying at least phase-encode 

and read gradient pulses following each RF pulse to acquire a 
set of data lines in k- space; 

a means (60) for changing a spectral offset with the 

generation of each data line sufficiently slowly that the 
steady-state condition does not need to be reestablished 
between the acquisition of consecutive data lines; 

a means (74) for combining steady-state signals 

derived from a plurality of data lines. 

28. The apparatus as set forth in claim 27, wherein 
the spectral offset is changed sufficiently slowly that the 
steady-state condition is. reversibly modulated between the 
acquisition of consecutive data lines. 

29. The apparatus as set forth in either one of 
claims 27 and 28, wherein a plurality of sets of data lines are 
collected at different spectral offsets and further including: 

a sorting means (70) for sorting the data lines of 

each data set to match up like data lines from each data set; 

a processor (78) for Fourier transforming the 

combined data lines into an image representation. 

30. The apparatus as set forth in any one of 
claims 2 7-29, further including: 

a means (36, 60) for adjusting an RF phase of a 

transmitter and a receiver at least piecewise continuously. 

3 1 . The apparatus as set forth in any one of claims 
27-3 0, wherein the RF pulses are applied in a pulse train, each 
RF pulse train containing at least three spectral offset 
components . 
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